Spin Conductance in one-dimensional Spin-Phonon systems 
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We present results for the spin conductance of the one dimensional spin- 1/2 Heisenberg and 
XY model coupled to phonons. We apply an approach based on the Stochastic Series Expansion 
(Quantum Monte Carlo) method to evaluate the conductance for a variety of phonon dispersions 
and values of spin-phonon coupling. From our numerical simulations and analytical arguments, we 
derive several scaling laws for the conductance. 
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Introduction — Controlling and understanding the dy- 
namics and transport of spins is an experimental and 
theoretical challenge aiming at the development of meso- 
scopic systems for spintronics applications. In particular, 
the interest is recently focusing on the effect of interac- 
tions and the coupling of spin systems to the environ- 
ment; several proposals for the generation and detection 
of spin currents exist i As a key quantity characteriz- 
ing ballistic mesoscopic spin transport is the (spin) con- 
ductance, it is important to understand its behavior in 
prototype interacting spin systems coupled to lattice vi- 
brations - phonons. 

Related to this problem, since the 60's, the spin- 
Peierls transition, the dimerization of a lattice due to the 
spin-phonon interaction in the bulk, has been discussed 
for quasi-one dimensional compounds as TTF, TCNQ, 
CuGeOa. A variety of analytical and numerical meth- 
ods are still employed in the study of the dependence of 
the transition on the phonon spectrum and spin-phonon 
coupling^ 

Finally, in close relation to this work, spurred by the 
experimental finding of an exceptionally high thermal 
conductivity^ in quasi-one dimensional transition metal 
compounds, there is recently considerable experimental 
and theoretical research interest in the finite temperature 
transport properties of spin models. 6,7,8 In particular, 
this activity leads to the development of novel numerical 
simulation techniques^ for the calculation of zero and fi- 
nite temperature transport properties in one dimensional 
electronic/spin models. 

In this paper we address the fundamental question of 
the spin conductance of the spin-1/2 Heisenberg model 
and the XY model coupled to phonons. Thus this study, 
on the one hand, sheds light from a new perspective on 
the old question of the spin-Peierls transition and on the 
other hand addresses a prototype question in mesoscopic 
physics; it is of direct experimental relevance, as a de- 
tailed experimental set-up for the observation of the spin 
conductance was recently proposed. 10 

Quantum Monte Carlo (QMC) methods have already 
been successfully applied to the evaluation of thermody- 
namic quantities in spin-phonon systems We will 
use here the Stochastic Series Expansion (SSE) which 
allows the calculation of the conductance of interacting 



spin systems^ in this original application we extend it to 
coupled spin-phonon systems. 

In the following we will first discuss the model and 
the method we are using, then we will present data for 
the conductance in the Heisenberg model for optical and 
acoustic phonon dispersions, and derive several scaling 
laws. Finally, we will show results for the XY model 
coupled to optical phonons. 

Model and method — The nominal spin-Peierls Hamil- 
tonian describing the coupling of the one dimensional 
Heisenberg model to lattice vibrations is given by 

H = [J + X(x n - X n +l)] (S n ■ S w+ i) 

+ ^o/2 [2p 2 n + (x n - x n+1 ) 2 ] , (1) 

n 

where S„ are spin-1/2 operators and x n , p n the displace- 
ment and momentum operators of the ions {u)q is the 
characteristic phonon frequency and here we set h = 1). 

Because of the minus sign problem, this system eludes 
a direct analysis using the QMC method that we intend 
to employ here. Hence, we resort to the following gener- 
alized Hamiltonian whic h is a n extension of the Hamil- 
tonian discussed in Refs. I."ll lL 
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H = 



^2 yJ + cti{a n + al) + a 2 (a n+ i + a] l+1 )j ■ h n 

[ W ( a n a n + ar\+l a n+l) - ^l(4an+l + 4+1°™) 



n=0 



w+ ( a n a l+i + a n +ia n ) 



(2) 



where h n = C + S^S^ +1 - (S+S n+1 + S n S , „ +1 )/2 and a n 
(ajj) are phonon annihilation (creation) operators. From 
Eq. one recovers the spin-Peierls model Eq. (JIJ if one 
sets 



A/V2: 



Oil 



LU 



2ui = 



(3) 



and C = 0, J = J. However, as a finite C < —1/4 is 
required for the implementation of the SSEi2iA4 we shift 
(following Ref. 0) the boson operators a n and J, 
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On a n -(C/2){a 1 +a 2 )/(LO-Lu 1 -LU + ) (4) 
J -» J+ (ai + a 2 ) 2 C7(w- wi - w+) (5) 

so that C* = C, J 7^ J . Still, to sample the Hamiltonian 
Eq. (J2J with the SSE we have to meet the requirement 
that all matrix elements of H are negative^*!* implying 
ai > 0, a 2 > 0, J > 0, or explicitly (C = -1/4): 

(ai + a 2 ) /& < uj — cjj — a-> + , ai,a2>0. (6) 

The conditions in Eq. © thus forbid a direct analysis of 
the spin-Peierls Hamiltonian Eq. Q as may be seen from 
Eq. J3J). However, as we will see we find two scaling laws 
which enable us to make predictions for the spin-Peierls 
model Eq. QJ. 

We define the operators j n and P nll 



3t. 



iJ x e 
2h 




a l [ 2co ] 



FIG. 1: The conductance as a function of ai for various uj < J 
for optical phonons (uji = = cj + and 0:2 =0). ai is given in 
units 2o». The dotted line is situated at the critical coupling 
of Ref. ITU The inset shows the same curves with a% in units 
of J. 



which are the current operator for magnetic moments at 
site n (we denote the "charge" quantum by e in analogy 
to the electric current) and the potential operator of a 
local "voltage drop" at site to, respectively. The spin 
conductance g is the linear response of the operator j n 
to the perturbation P m , 



q := lim Rc 



It is independent of n or m and will be evaluated accord- 
ing to Ref. while the SSE method for the Hamiltonian 
Eq. @ is implemented following Refs. IMTll Note that 
we are considering only low temperatures and the effect 
of leads on the conductance is not included. 

For the update of the spin degrees of freedom we use 
the very efficient scheme of "operator-loop updates"^ 
The simplest version of an update for the phonons is 
to suggest a specific change in the configuration and ac- 
cept/reject it with a certain probability. Here we want to 
emphasize that we obtain smaller autocorrelation times 
by using an appropriate cluster update scheme, reminis- 
cent of the loop update (see appendix, section ^ . 

In our simulations we use a system size N — 192, a 
temperature T = 0.02J/ks and typically perform 2 • 10 5 
MC sweeps. 

It should also be mentioned that in our simulations 
there is no need to introduce a bound for the phonon 
occupation number, but for technical reasons we impose 
a restriction of maximum 48 phonons per site; since it 
is much larger than the phonon numbers created during 
the simulation this introduces no further error. 

Optical Phonons — The two parameters uj\ and uj + 
lead to a nontrivial dispersion relation of the phonon 
modes. For the moment we will set them to zero and con- 
sider phonons of constant frequency uj. This (Einstein) 
model was discussed for ot\ — —a 2 using the DMRG 



methodic and for a 2 = with the SSE^i Both find 
that if a\ is larger than a critical coupling a c the sys- 
tem undergoes a dimerization transition. In Ref. the 
authors analyze different values for uj finding that a c is 
approximately linear in u> if u> < J. The authors in Ref. 
111! discuss only one phonon frequency, namely ui = J/8. 
They find that a c — l.3ui. (Data in our notation.) We 
will use this result since we focus on uj < J. The reason 
for this is that in many materials of actual interest J is 
significantly larger than the Debye temperature; we also 
find that our QMC method converges better for small uj. 

For the conductance g(T — 0) we expect the form of 
a step function: Above the critical coupling a% > a c it 
is zero due to the opening of the gap; for a\ < a c the 
system is supposed to be in a Luttinger phase, where 
the conductance is given by the Luttinger parameter K . 
Since the Hamiltonian Eq. is spin-rotationally invari- 
ant, the latter is fixed to the value of 1/2 (see Ref. IT^l . 

The fact that g should be constant for ol\ < a c is con- 
sistent with an effective Hamiltonian that we can derive 
using second order perturbation theory (cf. Refs. [TsHT^l 
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S n -S n +i — ^-j— S n -S n +2- (7) 



Indeed, for a 2 = we find no change in the conductance 
to order a\ since g does not depend on the bandwidth J. 

In Fig. ^ we present results for the spin conductance 
g for u)i = = uj + and a 2 — 0. Within error bars we 
confirm the expectation that the conductance remains at 
1/2 until the critical coupling is reached. In the gapped 
phase the conductance goes rather smoothly to zero in- 
stead of step-like which is due to the finite temperature 
T = 0.02J we are using. The point where the conduc- 
tance has decayed to zero gives a rough estimate for the 
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coupling aT where the gap satisfies A (cut) = ksT . (This 
can be tested using data from Ref. fl5L if we assume that 
A(ai, CK2) = A(ai — 02) — which is a natural assumption 
since this is the scaling relation of a c to be derived be- 
low, cf. Eq. (| 1 31) - For u> = J/20 we find ot=o.02J — o c ~ 
0.088 J, i.e., (o T =o.02j - o c )/o; a73 m 0.784 to be com- 
pared with Fig. [3 below.) 

Subsequent calculations we performed for T — 0.01J 
find an increased statistical error, but reveal no tem- 
perature dependence which implies that ot=o.ou an d 
«T=o.02jare very close to each other. Note that the gap 
has a typical Kosterlitz-Thouless form^ first it opens 
very slowly, but then it grows rather rapidly. 

Since g(T = 0) is a step function — which is fully de- 
scribed by the critical coupling where the conductance 
drops to zero — it inherits the scaling law of a c : 



g(ai, lu) = g(ai/u) for T = 0. 



(8) 
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FIG. 2: The conductance as a function of (qi — a. c )/uP' 7 ' i . 
The inset shows the same data as a function of ai — a c in 
units of J . 

To analyze the behavior of g in the gapped region more 
closely, it is instructive to plot g versus ct\ — a c . This is 
done in Fig. [21 We find good agreement with another 
scaling relation — different from Eq. (JSJ) — specific to the 
gapped phase, namely, 



7(0:1,0;) = g [(oi - a c )/u> K 



(9) 



where k = 0.73 ± 0.05. ft is natural to assume that in 
the gapped region g depends on the parameters of the 
system only through j3A. This would in turn imply that 
the scaling relation Eq. © holds for the gap A. Note 
that this relation is not consistent with the exponential 
behavior of A typical for a Kosterlitz-Thouless transition. 
However, we should expect Eq. © to be approximately 
valid for not too large a and uj. We note that the data for 
A(oj = 0.05,0.5) from Ref. ITEI show also good agreement 
with Eq. © when k = 0.6. 
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FIG. 3: The conductance as a function of ai (in units of 
Juj* /uj) (a.2 — 0, uj = J/8) for various dispersions u)(k). The 
dotted line indicates the critical coupling of Ref.fTll. The inset 
shows the same data when ai is given in units of J. 



Acoustical phonons — Now we will consider nontrivial 
phonon dispersions, wi / / lu + . Performing a canoni- 
cal Bogoliubov transformation in momentum (k-) space, 
we obtain the phonon frequency dispersion 



w(fc) = 2vV - lji cos(fc)] 2 - [oj+ cos(fc)] 2 



(10) 



We see that if u = u>i + co + , as, e.g. in the spin- 
Peierls model Eq. (JTJ, uj(k = 0) = 0, whereas if w > 
u>i + uj + there is a gap in the spectrum. This dispersion 
is illustrated in Fig. 




FIG. 4: The phonon dispersion Eq. IllOt for u>\ = w/2 and u> + 
is 0, (3/4)<^i, (15/16)o;i, u>i (from top to bottom). At u> + — 
the parameter oij gives the width of the dispersion and uj its 
mean value. 

After diagonalizing the phonon part of the Hamilto- 
nian Eq. © we proceed similarly with the spin-phonon 
coupling term. To this end, adopting for convenience a 
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fermion representation (S + i— > cjj by the Jordan and 
Wigner transformation, the phonon-fermion scattering 
term 

n 

is written in k-space as 

Since we are here interested only in the low temper- 
ature regime we may assume that the variables k, q lie 
at the Fermi momenta ±7r/2. The process responsible 
for the decrease of the conductance should then be the 
scattering between two Fermi points by a phonon with 
momentum 7r and the only relevant frequency should 
be W (pi). Thus, for a general dispersion relation, we 
should recover the same result for g as for the opti- 
cal phonons of the prev iou s section, by m apping u — > 
u)(n)/2, aii — > aiy/id(n)/ y/2(id + w\ + lu + ). We expect 
then the conductance g to depend only on a\/iu* with 
lu* = y/(io + id 1 + uu + )uu(tt); to verify this expectation we 
show g in Fig. El for a set of parameters u>, oj\, lu + com- 
patible with Eq. ©. 

The figure indeed suggests that g(ai,id,u>i,u) + ) = 
g(otiid/(d* and together with Eq. JSJ we arrive at 
the desired scaling relation 



g(ai, lu, u>i, lu + ) — g(V2ai/uj*), for T = 0. (11) 




FIG. 5: The conductance as a function of ai — Q2 for various 
0.2/0.1. (The inset shows g vs ct\ + 02- The order of the curve 
-top to bottom- is the same as in the legend.) Here lu = J/8, 
wi = = lu + , T = 0.02J/fc s . 

Scaling with ct\ and a.2 — Let us again consider the 
Hamiltonian in k-space as in the previous section. If we 
have <x\ and cti different from zero, the coupling terms 
read 

\ £ A k , q ( ai + e^a 2 )c\c q (a^ q + a\_ k ). (12) 

k,q 



Here Ak, q is a prefactor — whose precise form is for the 
following discussion irrelevant. At low temperatures we 
have again k,q = ±n/2 which implies that k — q = 0, tt. 
Hence only two phonon modes are important: u>(0) and 
u>(ii). The w(0) mode describes scattering on one Fermi 
point which means that it probably rescales the parame- 
ters of the system, but has otherwise no influence on the 
transition. For the phase transition only the ui(ir) mode 
is important. As we see from Eq. I|12|) the mode with 
ui(0) couples to the spin degrees of freedom via a± + a 2 
and the mode u>(tt) via a\ — a 2 . As the decay of the 
conductance is due to the opening of the gap, we are led 
to assume that g depends on ol\ — a 2 , but not on ai+a 2 , 
i.e., 

g(ax,a.2) = g(a\ — 012)- for T = (13) 

(Alternatively, one can motivate this scaling law with 
a Mean Field Theory ansatz which is outlined in the ap- 
pendix, section O) 

In Fig. we confirm this scaling by the overlap of 
curves as a function of ct\ — a 2l while in the inset we 
indeed find that the dimerization transition is delayed 
by increasing a 2 /a\ 1 till a 2 — ct\ where it is absent as 
predicted by Eq. 1(15)1 . 

A further test of Eq. I|13|) is provided by comparing 
data for a c . For lu = J/8 and a 2 = Ref. [Il] finds 
a c = 0.163J (with 4% uncertainty), whereas according 
to Ref. for the same lu but u\ = — a 2 the transition 
occurs when ol\ — a 2 = 0.176 J. The difference of the two 
values is 8% which means reasonable agreement. 
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FIG. 6: The conductance of the XY model as a function of 
oti {a 2 = 0) for two u (here T = 0.01J). The dotted line 
gives the critical coupling if one assumes that the scaling law 
a c (ai,ot2) ~ a c (ai — a 2 ) holds. The dot-dashed line gives 
the critical a c extracted from Fig. [H] 

XY model — In this section we turn to the XY model, 
that is we replace S„ • S n +i in the Hamiltonian Eq. by 
S%Sn +1 + S*S^ +1 . The model was studied by Caron and 
Moukouri 18 for ct\ — —a 2 . They found a phase tran- 
sition similar to the one in the spin-Peierls model but 
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FIG. 7: The conductance of the XY model as a function of 
ai-Q2 (solid curves), (qi — Q2) 1,2 /( a i +Q2) ' 2 (dashed curves 
which are offset by —0.07 for clarity) and ot\ + 0:2 (inset) for 
various 012/011 with ui — J/8, T — 0.01J. 
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We computed the conductance for two 



phonon frequencies. The results are shown in Fig. [HJ We 
find that the zero temperature conductance obeys the 
scaling law 



g{ax,u) = g{a x /u ■ ) 



(14) 



It is important to point out that the XY model is not 
spin-rotationally invariant and therefore the value of the 
conductance is not fixed in the Luttinger phase, but 
rather decreases monotonously, as the Fig. shows. In 
the gapped phase the conductance decays again quite 
smoothly to zero, such that it is barely possible to decide 
where exactly the transition takes place. 

We therefore employed here two different methods for 
calculating a c . First QMC simulations for the staggered 
kinetic energy which give a rigorous upper bound (details 
are deferred to the appendix, section lB|l . second we used 
the data for a c from Ref. [l8| while assuming that the 
scaling law Eq. (|13fl . g(a±, a-i) — g{a\ — a 2 ), holds 

The two results for a c are indicated by lines in Fig. [B] 
for comparison. The difference between the two predic- 
tions for a c is about 15%. One should emphasize again 
that both methods are not numerically exact: the QMC 
yields only an upper bound for a c , and the other method 
relies on the validity of Eq. I|13[) which is examined in 
Fig. [7] by considering nonzero a 2 . For the parameter ac- 
cessible to us Eq. (|13f) seems to be approximately, but 
not exactly fulfilled. 

In the Heisenberg model the dependence on ax, a 2 and 
li was described by two separate scaling laws Eqs. ijSjl 
and (|13fl . For the XY model the situation is more com- 
plicated, as we found that for a 2 = «i/4 the scaling 
law Eq. (|14fl is violated. This implies that a general- 
ization of Eq. I|13|) which holds for the XY model must 
depend on li. We will not try to find such a general- 
ization here, but want to point out that for the special 



case of lu = J/8 one can obtain a much better fit for our 
data (for a^/ati = 0, 0.25, 0.5) than Eq. (|T3|) by using the 
following modified scaling ansatz (see Fig. 0) 



g(ax,a 2 ) = g 



0L2) 



Oi\ — a 2 

Oil + Q'2 



Conclusions — We have shown that the known scal- 
ing laws for a c in the Heisenberg and the XY model, 
i.e., Eqs. (JSJ & (|14fl also hold for the conductance. For 
the spin-rotationally invariant case we found a general- 
ization for acoustic phonons Eq. I|ll(l as well as a scaling 
relation for the gap Eq. © . We also discussed the differ- 
ence between two possible couplings to phonons (a 2 = 
& a 2 = —oil). For the spin-rotationally invariant case 
the two ways of coupling are related by a simple scal- 
ing law Eq. (|13fl . whereas for the XY chain the situa- 
tion appears to be more complicated. Assuming that our 
scaling laws are not restricted we may make a predic- 
tion for the conductance of the spin-Peierls model Eq. 
(JTJ by putting together Eqs. ©, CU an d D31 as weu 
as data for the critical coupling from Ref. [IJ We ex- 
pect a step-like conductance which drops to zero when 
A = 0.65- 2 5 / 4 w ~ 1.55w - 

One final comment should be made. Over the last few 
years, it has generally agreed upon that the conductance 
measured in experiments corresponds to the conductance 
of a system coupled to noninteracting leads. The leads - 
which were not considered here - may strongly affect the 
conductance>i2i2S We can argue however that the dimer- 
ization transition we observed will drive the conductance 
to zero even if it occurs only locally, i.e. in the interact- 
ing part system and not in the leads. Hence, we think 
that our results have a clear experimental consequence. 

We would like to acknowledge useful discussions with 
D. Baeriswyl and financial support by the E.U. grant 
MIRG-CT-2004-510543. 

APPENDIX A: MC UPDATE FOR THE 
PHONONS 

We write the Hamiltonian Eq. in the form H = 

hn^ ,k \ where n is the site (bond) index and indices 
(i,j,k) denote different types of (local) operators. 

The parameter j (resp. k) assumes values from — 1 
to 1 and refers to the change of the phonon number on 
site n (resp. n + 1) resulting from the action of the op- 
erator hn^'^. In the sequel we will refer to j or k as 
phonon parameters. The index i refers to the action on 
the spin degrees of freedom. For example, the operators 
for different i and j = = k read 

= JS+S- +1 /2 h^ = (h^ s >) 1 
h (3,o,o) = C +JS*S* +1 +Li(ata n + al +1 a n+1 ). 
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Following Ref. 0, we expand the Boltzmann weight 
in terms of the local operators, i.e., we sample the op- 
erator strings nPCIMOR + ^) according to their weight 
in the Boltzmann factor are states in the S z -phonon 
occupation-number basis and h(l) is the local operator 

hn^^ which appears on the Zth position of the string). 
The spin degrees of freedom can be updated using the 
very efficient scheme of operator loops^ During this up- 
date the parameters j and k remain unchanged. All that 
remains to be done is to devise an additional update for 
the phonon parameters. 

This is achieved as follows: As the operators a n and a J, 
appear always pairwise in the string, we have to increase 
the phonon parameter at one position (of the string) and 
decrease it at another position. In principle the two posi- 
tions can be chosen at random, but here we want to point 
out that this can be done more efficiently: We choose an 
operator hn in the string. Next, we decrease (resp. 
increase, depending on an appropriate probability) one 
of its phonon parameters (associated with a site which 
we will call m). For the previous (resp. next) operator in 
the string we choose — with a certain probability — one of 
three options: either (i) we update the phonon number 
on site m by +1 or (ii) we cancel the operation and re- 
store the initial configuration, or (iii) we increase (resp. 
decrease) the phonon parameter related with site m. In 
case (ii) or (iii) the procedure stops, while in the first 
case it is iterated until we finally choose one of the other 
options. 

Even though one might suspect that because of option 
(ii) the algorithm is inefficient — since it is reminiscent 
of the bounce process in Ref. Q — we find that our new 
update scheme reduces the autocorrelation time signif- 
icantly. The reason for the latter may be seen in the 
fact that we choose only one position of the two phonon 
operators at random, the other is — in some optimized 
way — chosen by the MC algorithm. 
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FIG. 8: Absolute value of the contributions of the staggered 
kinetic energy (|ftics|) (see text) for the XY-model with u = 
J/8 at T = 0.01J. Here we use periodic boundary conditions 
and I0 4 MC sweeps. 
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FIG. 9: Absolute value of the contributions of the staggered 
kinetic energy (l^ryl) (see text) for the Heisenberg-model 
with lo = J/8 at T = 0.01J. The dashed line give the crit- 
ical coupling from Ref. ITU Here we use periodic boundary 
conditions and 10 4 MC sweeps. 



APPENDIX B: EVALUATION OF THE 
CRITICAL COUPLING WITH QMC 

The simplest method to determine a c with Monte 
Carlo is measuring the order parameter. In our case we 
have two options for the order parameter: one is 

n 

the staggered bond lengthj2I another is the staggered ki- 
netic energy^ 

**y = £(-!)"(# S n+i + S-S+ +1 )/(2N). 

n 

Since these operators are linear combinations of terms 
appearing in the Hamiltonian, there is a simple method 
to evaluate the expectation values with the SSEii^ To 
this end we define n a>n (and n xv ,n) which are the sum 
of the numbers of appearances of the operators Ca n and 
Cat (J/2S+S- +1 and J/2S-S+ +1 ) in the string of opera- 
tors created and updated in the simulation process. The 
expectation values of the order parameters are related 
with the mean values of the numbers n 0j „ and n xy>n ; the 
relation being 

(*«) = £(-i) n <-»«,»>/(JW 

n 

and 

(**v> = 52(-l) n (-n*y,n)/(NpJ). 

n 

Since the contributions to the order parameters may 
be positive or negative, it is difficult to detect a finite 
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dimerization — especially, when we are close to the criti- 
cal coupling. Hence, it is difficult to decide whether we 
are in the gapped phase or not. We therefore calculate 
the following expectation values 

(i*«i> = (|E(- 1 ) nn ».» 



and 



(One should note that in general (|$ a |) 7^ (|$a|) where 
l^al = V^a-) These have non-zero values, but in the 
gapless phase they should go to zero in the thermody- 
namic limit. We note that while ($ a ) is independent of 
C the expectation value (|$ a |) does depend on the shift 
C. Because of this and the fact that (|$ a |) > (l^xyl) 
(which suggests that (|$a:y|) decays faster with N) we 
prefer the second order parameter to the first and will 
only discuss (|$aT/|)- 

Results for the staggered kinetic energy, are 

presented in Figs. (El and E| For the XY model (Fig. EJ 
the order parameter (l^xj/l) decreases with system size as 
long as a>i < 0.475J. For a% > 0.5 J it seems to increase 
with system size rather than decrease. We conclude that 
a c < 0.5 J. A correct assessment of the critical coupling 
remains difficult. A rough estimate for a c is obtained by 
assuming that a c lies in the interval 0.475 J < a± < 0.5 J 
where the behavior of (l^xj/l) changes. This yields a c « 
0.4875 J ± 0.0125 J. In principle, this estimate gives only 
an upper bound for a c , but when we apply the same 
method to the spin-rotationally invariant case, we find 



good agreement with the value for a c from Ref. [Tl| (see 
Fig.EJ. 



APPENDIX C: MEAN FIELD 

To better understand the physics of the model a mean 
field analysis is helpful. The main idea is to replace the 
displacements (a n + aj,)/\/2 by classical variables x n . A 
common but general ansatz for these variables is x n — x+ 
(— l) n A. Substituting this ansatz into the Hamiltonian 
Eq. J2J) we obtain, 



Hmf — E JnSn ■ S n+ i (CI) 

n 

J n = J + V2(ai + a 2 )x + (-l) n y/2{a l - a 2 )A. 

One sees that the spin-Peierls model Eq. (JIJ does not 
depend on x because of a\ = —a 2 . At the mean field 
theory level we arrive at two conclusions, (i) there is a 
shift in J proportional to (a% + a 2 ) and (ii) there is a 
dimerization transition where the dimerization is propor- 
tional to (cki — 02). In the spin-Peierls model (ai = — a 2 ) 
the shift in J is absent, whereas in a toy model with 
Oil = a 2 no dimerization should be found. In our simu- 
lations (ai, ot2 > 0) we always have a shift but, since the 
conductance at T = does not depend on J, we expect 
that g is independent of a\ + a 2 and we arrive at the 
scaling law 



g(ai,a 2 ) = g(ai - a 2 ). 



(C2) 
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